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ABSTRACT 

GRB 090709A is a long gamma-ray burst (GRB) discovered by Swift, featuring a bright X- 
ray afterglow as well as a faint infrared transient with very red and peculiar colors. The burst 
attracted a large interest because of a possible quasi-periodicity at P = 8.1 s in the prompt 
emission, suggesting that it could have a different origin with respect to standard, long GRBs. 
In order to understand the nature of this burst, we obtained a target of opportunity observation 
with XMM-Newton. X-ray spectroscopy, based on XMM-Newton and Swift data, allowed us 
to model the significant excess in photoelectric absorption with respect to the Galactic value 
as due to a large column density (~6.5 x 10 22 cm -2 ) in the GRB host, located at z ~ 4.2. 
Such a picture is also consistent with the infrared transient's properties. Re-analysis of the 
prompt emission, based on International Gamma-Ray Astrophysics Laboratory and on Swift 
data, excludes any significant modulation at P = 8.1 s. Thus, we conclude that GRB 090709A 
is a distant, standard, long GRB. 
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The bright, long gamma-ray burst (GRB) GRB 090709 A was dis- 
covered by the Burst Alert Telesc ope (BAT) onboard Swift on 2009 
July 9 at T = 07:38:34.59 UT dMorris et alj|2009h . The prompt 
emission had a complex structure (see Fig. [I), with a multi-peaked 
light curve lasting tgo ~ 89 s. The peak flux was 7. 8±0. 3 ph cm -2 
s" 1 at To + 21 s and the fluence w as ( 2.57±0.03) x 10"^ erg cm" 2 
in the 15-350 keV energy range dSakamoto etai]|2009l) . A bright 
X-ray afterglow was observed by the X-ray Telescop e (XRT) on- 
board Swift starting as soon as 77 s after the trigger (Mor ris et al.l 
2009). Analysis of XRT data up to To + 0.3 days yielded evidence 
for a break in the decay at To + 0.1 days. The XRT spectrum was 
described by a power law with an absorbing column exceeding by 
a factor ~3 the Galactic va lue in the direction of GRB 090709A 
dRowlinson & Morrisll2009l) . 

On the optical/infrared side, several follow-up observations 
of the field of GRB 090709A were performed. A possible, faint 
transient with very red colors was detected in early observa- 
tions in the near infrared, performed within a few minutes from 
To by automated instruments such as the PAIRITEL Telescope 
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(Moraan, Bloom & Kleii] [ 2009 ) , the Palomar Observatory's 60- 
inch telescope (Cenkoetal. 2009) as well as the Faulkes North 
Telescope ( Guidorzi e t al.l2009h . The peculiar colors of such a tran- 
sient prompted to suggest a very high redshift (z ~ 10) for such 
GRB. The same infrared source was possibly d etected in a Sub - 
aru image collected ~2 hours after the trigger dAoki et al]|2009h . 
A very deep observation performed with the 10.4 m Gran Telesco- 
pio CANARIAS telescope ~41 hours after the GRB detected no 
P 1 

dCastro-Tirado et alj|2009l) 

GRB 090709A attracte d much interest because o f a very pecu- 
liar timing phenomenology, llvlarkwardt et alj d2009l) reported evi- 
dence for quasi-periodical pulsations at ~8 s in the prompt emis- 
sion, based on analysis of Swift/BAT data in the 15-350 keV energy 
range. Such a quasi-periodicity was the n confirmed by the anal - 
ysis of Konus-Wind and Konus-RF data dGolenetskii et alj 2009). 
as well as by analysis of INTEGRAL/ACS data dG6tz et alj|2009f) . 
Taken at face value, such a result is tantalizing. This would be the 
first detection of periodic activity in the prompt gamma-ray emis- 
sion of a GRB, which could yield rare information on the behavior 
of the central engine. 

Indeed, the pecul iar quasi-periodicity im medi ately prompted 
sever al authors (e.g. Ma rkwardt et al.l 120091 and iGuidorzi et al.l 
120091) to suggest that GRB 090709A could be different in origin 
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Figure 1. Light curve of GRB 090709A obtained with the Swift/BAT in the 
15-350 keV. The time binning is 0.5 s. 



with respect to standard long GRBs, and possibly related to activ- 
ity of a magnetar, either within the Galaxy, or extragalactic. How- 
ever, a search for pulsations i n early afterglow data col lected by 
Swift/XKT yielded null results dMirabal & Gotthelj|2009l) . 

To shed light on the nature of the peculiar GRB 090709A we 
asked for an XMM-Newton Target of Opportunity (ToO) observa- 
tion, aimed at characterizing the afterglow emission at a later stage 
(about 2 days after the event), when the target is too faint to be stud- 
ied with Swift/XKT. Here we report on the outcome of the XMM- 
Newton ToO, which, coupled to a re-analysis of Swift and INTE- 
GRAL data, allowed us to obtain a comprehensive description of 
the high-energy phenomenology of GRB 090709A. 



2 IS THE X-RAY AFTERGLOW "PECULIAR"? 

2.1 The XMM-Newton view of GRB 090709 A 

The ToO observation started on 2009 July 11 at 07:56:50 UTC 
(48.3 hr after the burst) and lasted 24.5 ks. All the European Pho- 
ton Imaging Camera (EPIC) detectors were operated in Full Frame 
mode (imaging across the whole field of view, with a time resolu- 
tion of 73 ms and 2.6 s in the pn and in the two MOS cameras, 
respectively), using the thin optical filter. We processed Observa- 
tion Data Files using the most recent release of the XMM-Newton 
Science Analysis Software (SASv9.0). 

No significant particle background episodes affected the ob- 
servation. The afterglow is clearly detected in all of the EPIC cam- 
eras. The position, RA = 19 h 19 m 42!6, Dec. = +60°43'35'.'0 
(J200 0), with a la error of 1'.'5, is consistent with the XRT posi- 
tion dOsborne et"ai]|2009h . 

Source photons were selected from a circular region (30 arcsec 
radius) centered on the target. Background events were extracted 
from a source-free region in the same chip as the target. The source 
average background-subtracted count rate in the 0.3-8 keV energy 
range is 0.155 ± 0.003 counts s"\ 0.051 ± 0.002 counts s* 1 and 
0.048 ± 0.002 counts s" 1 in the pn, MOS1 and MOS2 cameras, 
respectively. Background contributes ~5% more counts in the ex- 
traction region. 

2. 1. 1 Search for pulsations 

As a first step, after correcting photon arrival times to the solar 
system's barycenter using the task BARYCEN, we searched for pul- 
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Figure 2. Results of the timing analysis on EPIC/pn data (0.5-10 keV), 
collected ~48 hours after the trigger. Upper panel: the power spectrum is 
shown together with the threshold for the detection of sinusoidal signals at 
the 3er confidence level. Lower panel: upper limits on the pulsed fraction. 
See text for details. 

sations in th e GRB afterglow usi ng two methods. First, we used 
the Z 2 test teuccheri et alj[l983l) . with the number of harmonics 
n being varied from 1 to 2. Second, a Fourier analy sis of the light- 
curves was performed using the method described in Israel & Stella 
(1996). The analysed period range spans from 150 ms up to 10 s 
(~262 000 total period trials). No significant (periodic or quasi- 
periodic) signal was found with either method, searching the whole 
period range or restricting the search around the 8 s hypothetical pe- 
ri od. The upper l i mits o n the pulsed fraction, computed according 
to lVaughanetalJ l ll994l) . are shown in Figure 2 as a function of the 
frequency (the curves refer to the 3<r upper limits on the sinusoid 
semi-amplitude pulsed fraction in the 0.5-10 keV energy range). 
Upper limits at the 8 s hypothetical period range between 26% 
when considering all the Fourier frequencies in the whole spec- 
trum and 23% for the narrower search (see insets in Figure[2] lower 
panel). 

2.1.2 Spectral analysis 

Source spectra were rebinned in order to oversample the instrumen- 
tal energetic resolution by a maximum factor 3, or to have at least 
25 counts per bin. Ad-hoc response matrices and effective area files 
were generated using the SAS tasks RMFGEN and ARFGEN, respec- 
tively. 

Spectral analysis was performed using the XSPEC software 
(vl2.4.0). Errors on spectral parameters are given at the 90% con- 
fidence level for a single parameter of interest. Simultaneous mod- 
elling of pn, MOS1 and MOS2 spectra with an absorbed power 
law model results in a rather poor fit {xl = 1-52, 152 d.o.f.). 
The best fit photon index is V = 2.15 ± 0.05, while the absorb- 
ing column (JV H = (2.5 ± 0.2) x 10 21 cm -2 ) is significantly 
larger than the Galactic value in the direction of GR B 090709A 
(~6.7 x 10 20 cm -2 , according to lDickev & Lockmanll99d) . Such 
result are in broad agreem ent with the XRT ones reported by 
iRowlinson & Morris! J20091) . 

We then added to the Galactic absorption component (with 
JVh fixed to 6.7 x 10 20 cm -2 ) a redshifted absorption component 
Nh,z (abundances were set to Solar values). This model yields a 
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Figure 3. Soft X-ray light-curve of GRB 090709 Aspanning the time range 
from T + 77 s to To + 20 days. Blue triangles: Swift/XRT data; Red 
circle: XMM-Newton data. The broken power law best describing the flux 
decay is superimposed. Time intervals for XRT time-resolved spectroscopy 
are marked as A, B, and C. 

much better description of the data (xt = 1.11, 151 d.o.f.). The 
resulting photon index is T = 2.00 ± 0.05, The fit gives a very 
large intrinsic column density TVh.z = (9.7 ± 2.0) x 10 22 cm~ 2 , 
while the redshift is constrained to 5.1 ± 0.4. The time-averaged, 
observed flux is 5.6 x 10 -13 erg cm' 2 s _1 in the 0.5-10 keV en- 
ergy range. 

2.2 Swift/XRT data 

In order to put the XMM-Newton observation in context, we 
retrieved and analysed the XRT observations of the afterglow. The 
XRT CCD detector (0.2-10 keV) started observing the field of 
GRB 090709A only 77 s after the BAT trigger. Table □ reports 
the log of this and the subsequent XRT follow-up observations, 
performed in both photon counting (pc) and windowed timing (wt) 
modesQ that were used for this work (some wt observations last 
only a few seconds and we did not include them in our analysis). 

The data were processed with standard procedures 
(XRTPIPELINE version 0.12.3), filtering, and screening crite- 
ria by using FTOOLS in the HEASOFT package (ver. 6.6). For the 
timing and spectral analyses, we extracted the pc source events 
from a circle with a radius of 20 pixels (one pixel corresponds 
to about 2'.' 36) and the wt data from a 40 x 40 pixels box along 
the image strip. To estimate the background, we extracted pc 
and wt events from source-free regions far from the position of 
GRB090709A. 

The overall light curve of the afterglow of GRB 090709A 
is shown in Fig. |3] Count rates (for both 5w;/f/XRT and XMM- 
NewtonfEPlC) have been converted into fluxes using the best-fitting 
spectral models described below and in the previous section. The 
flux decays as a broken power law. The break occurs at 0.26 ± 0.05 
days after the GRB, when the index of the decay changes from 
-1.15 ±0.01 to -1.48 ±0.05. 



1 In pc mode the entire CCD is read every 2.507 s, while in wt mode only 
the central 200 columns are read and only one-dimensional imaging is pre- 
served, achieving a time resolution of 1.766 ms (see iHill et alj|2004l for 
more details). 
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Figure 4. Same as Fig. [2] using the data collected with Swift/XRT (0.5-10 
keV) in WT mode in the early afterglow phase (from To + 77 s to To + 539 
s). 

2.2.1 Search for pulsations 

We calculated the power spectrum for the first Swift dataset during 
which the time resolution was higher (WT mode) and with large 
enough statistics. The analysed period range spans from 4 ms up 
to 10 3 s (~262 000 total period trial) approximatively. Also in this 
case no significant signal was found searching the whole period 
range or restricting the search around the 8 s signal. Due to the 
presence of low-frequency noise (introduced by the source rapid 
decay during the Swift observation), meaningful upper limits in the 
range 60%-50% have been inferred only for the narrow search. The 
subtraction of high-order de-trending polynomials (taking account 
for the source rapid decay) does not change significantely the above 
results. We note the presence of a low-significance (~2rj) QPO-like 
feature in the power spectrum around lis (Fig. 4). 

2.2.2 Spectral analysis 

We extracted time-resolved spectra of the afterglow. In order to 
have a good photon statistics, we combined data collected within 
the three time intervals marked as A, B, and C in Fig. [3] 

A simple, absorbed power-law model does not reproduce well 
the spectra and results in a Nn column significantly larger than the 
Galactic value. As in the XMM-Newton case, we added a redshifted 
absorption component, which yielded a much better fit to the data. 
Results are reported in Table|2] 

In order to assess the possible spectral evolution as a function 
of time, we generated confidence ellipses for z vs. Nu,z and for 
F vs. Nh,z for the three datasets having the largest photon statis- 
tics (XRT A, XRT B and EPIC). As shown in Fig. [5] the strong 
correlation between the Nn,z and z parameters prevents from con- 
cluding that we are observing different values in different epochs. 
Conversely, a significant softening of the power law component is 
apparent. The photon index varies from ~ 1.7 in the early afterglow 
phase to ~2 at the time of the XMM-Newton observation. 

Indeed, fitting the above model simultaneously to all XRT and 
EPIC spectra, leaving the power law as the only spectral component 
free to vary as a function of the epoch, yields a very good result 
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Table 1. Journal of the Swift/XKI observations. The Swift/BAJ trigger time is 2009-07-09 07:38:35 UT. 



Sequence/Mode 


Start/End time (UT) 


Exposure" 


Net average count rate b 




yyyy-mm- 


dd hh:mm:ss 


(ks) 


(counts s~ ) 


00356890000/wt 


2009-07-09 07:39:52 


2009-07-09 07:47:33 


0.4 


35.4 ±0.3 


00356890000/pc 


2009-07-09 08:47:09 


2009-07-09 15:14:03 


8.9 


0.744 ± 0.009 


00356890001/pc 


2009-07-09 16:56:40 


2009-07-09 02:29:10 


7.9 


0.103 ± 0.004 


00356890002/pc 


2009-07-10 02:34:55 


2009-07-10 20:40:56 


10.2 


0.036 ± 0.002 


00356890259/wt 


2009-07-11 04:14:24 


2009-07-11 15:39:27 


12.0 


0.008 ± 0.003 


00356890003-4-5/pc 


2009-07-19 06:56:06 


2009-07-21 23:38:57 


30.4 


0.0011 ± 0.0002 


00356890006-7-9- 10-1 1/pc 


2009-07-22 02:39:01 


2009-07-28 22:45:58 


53.4 


0.0012 ± 0.0002 



a The exposure time is usually spread over several snapshots (single continuous pointings at the target) during 

each observation. 

b In the 0.5-10 keV energy band. 



Table 2. Results of time-resolved spectroscopy with Swift/XRT and XMM-Newton/EPIC. 





Swift/XKI A 


Swift/XSJ B 


Swift/XKT C 


XMM-Newton/EPIC 


Exp. Time (ks) 


0.4 


8.9 


17.9 


18.5 a /22.6 6 


Mean epoch (MJD) 


55021.322 


55021.500 


55022.248 


55023.480 


AT H (cm -2 , fixed) 


6.7 x 10 20 


6.7 x 10 20 


6.7 x 10 20 


6.7 x 10 20 


N H , Z (cm" 2 ) 


(6.2 ± 1.1) x 10 22 


(8.7 ± 1.7) x 10 22 


(7.lit*) x 10 22 


(9.7 ±2.2) x 10 22 


z 


4.1 ±0.3 


4.3 ±0.3 


3.9 ± 1.0 


5.1 ±0.5 


r 


1.69 ±0.03 


1.84 ±0.05 


2.05 ± 0.14 


2.00 ± 0.06 




1.13 


0.98 


1.12 


1.12 


d.o.f. 


291 


172 


39 


151 


F c (erg cm -2 s — 1 ) 


1.6 x 10" 9 


3.7 x 10" 11 


1.6 x 10" 12 


5.6 x 10" 13 



l b pn/MOS. 



photons' time of arrival to the solar system barycenter, we gener- 
ated a background-subtracted light curve in t he 15-150 keV energ y 
range using the mask-weighting technique dSenziani et al.|[2007l) . 
Then, we used a procedure similar to that already discussed for 
XMM-Newton and Swift/XKT data. We searched for possible sig- 
nals in the power spectrum. In this case we considered the whole 
set of Fourier frequencies. We searched for significant signals in 
the range from 20 ms to 50 s, approximately. We also included 
high-order de-trending polynomials to "rectify" (trend subtraction 
mode was selected) the light curve and minimize the low-frequency 
noise. Fig. [6] shows our best-case power spectrum, obtained after 
de-trending the light curve by using a third order polynomial. No 
significant signal was found, while a 3a upper limits in the 15%- 
25% range was obtained in the period interval between 5 s and 15 
s. We note that a peak at about 0.125 Hz, consistent to the reported 
period of 8 s, is present in the spectrum. The corresponding power 
estimate is very close to the 3<r detection threshold. 

We also studied the time-averaged spectrum of the prompt 
emission, generating a bac kground-subtracted s pectrum with the 
mask-weighting technique ( Senziani et al. 2007), considering the 
time interval from Tq to To + 150 s. The best fit model is a power 
law with photon index T = 1.28±0.02. The fluence in the 15-150 
keV energy range is 2.3 8 x 10 erg cm" . O ur results are fully 
consistent to the ones bv lSakamoto et al] d2009t) . 



c Observed flux in the 0.5-10 keV energy band. 

(Xt = 1-08, 671 d.o.f). Such an exercise yields a best fit redshift 
z = 4.2 ± 0.2 and intrinsic absorption AT HjZ = (6.5 ± 1.5) x 10 22 
cm -2 , while the epoch-dependent values of F and of the power- 
law normalization are very similar to the ones reported in Table[2] 
Taking into account the correlation between z and Nk,z, as well as 
the dependence of such parameters on the value of the Milky Way 
column densit)0, the 90% confidence level interval turn out to be 
3.7 < z < 4.5 and 4.8 x 10 22 cm" 2 < N H ,z < 7.7 x 10 22 cm -2 . 
We also note that such results are computed assuming Solar abun- 
dances in the redshifted absorber model. 



3 ANY PERIODICITY IN THE PROMPT EMISSION? 

Facing with a rather standard X-ray afterglow, we have studied the 
prompt emission in order to assess the significance of the claimed 
quasi-periodical variations. 



3.1 SwiWBAT data 

The coded-mask Burst Alert Telescope (BAT) onboard Swift col- 
lected event data (15-350 keV, with 0.1 ms time resolution) from 
GRB090709A in the time range from T - 170 s to T + 612 s. 
We retrieved the event list from the Swift archive. After correcting 



2 IStrattaetalJj2004l) showed that a reasonable guess of the uncertainty af- 
fecting Galactic column density estimates is ~ 30% at ~ 90% confidence 
level. 



3.2 INTEGRAL SPI/ACS data 



The Anti-Coincidence Shield (ACS) of the Spectrometer on INTE- 
GRAL (SPI) is routinely used as a nearly omni-directional detector 
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Figure 5. Upper panel: confidence contours (solid: 68%, dashed: 90%, dot- 
ted: 99%) for the redshift z vs. the intrinsic column density iVj] g. Lower 
panel: confidence contours (solid: 68%, dashed: 90%, dotted: 99%) for the 
photon index V vs. the intrinsic column density Nn z . 



for gamma-ray bursts dvon Kienlin et alj|2003l) . besides serving its 
main function as a veto for the SPI spectrometer. The ACS provides 
light curves binned at 50 ms, but without energy and directional in- 
formation. The low energy threshold is about 80 keV. 

GRB 090709A was located at an angle 9 = 70° from the SPI 
pointing direction, resulting in an optimal response for the ACS, 
which is most sensitive for directions orthogonal to the satellite 
pointing axis. 

In order to look for a possible periodic signal in the ACS 
data we perfomed the same analysis as in the Swift/BAT case (see 
above). The results are shown in Fig. [7] also in this case the search 
was negative. The 3a upper limits to any pulsed signal around 
0.125 Hz is 20%, and 6% and 8% for the blind search and the nar- 
row search in the 5-15 s period interval, respectively. 



4 DISCUSSION AND CONCLUSIONS 

4.1 The afterglow of a distant, standard, long GRB 

Prompted by early reports of the detection of a strong quasi- 
periodical signal in the prompt emission of GRB 090709 A, we have 
performed a ToO observation with XMM-Newton. 
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Figure 6. Similar to Fig.f2]for the prompt emission, using SwiftfBPCI data 
(15-150 keV) in the time range from To s to To + 100 s. The uppermost 
line (stepped) marks the threshold for the detection of sinusoidal signals at 
the 3<r c.l., while the lower one shows the corresponding upper limits on the 
pulsation amplitude. 
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Figure 7. Similar to Fig.|6]but for the prompt emission as seen by INTE- 
GRAL/SPl- ACS. The uppermost line (stepped) marks the threshold for the 
detection of sinusoidal signals at the 3cr c.l., while the lower one shows the 
corresponding upper limits on the pulsation amplitude. 



The afterglow of GRB 090709A as seen in soft X-rays with 
the EPIC instrument, is fully similar to other cases of typical, well- 
behaving, long GRBs. No pulsations are seen in the X-ray emis- 
sion ~48 hr after the trigger, with a 3a upper limit of 15-25% 
on the pulsed fraction, assuming a sinusoidal pulse shape. The en- 
ergy spectrum has a typical power-law shape, with a remarkable 
excess in photoelectric absorption (by a factor ~4) with respect to 
the Galactic JV H in the direction of GRB 090709A. 

Such a picture of a "standard" X-ray afterglow is completed 
by Swift/XKF observations, which started as soon as 77 s after the 
burst, and extended up to 20 days after the trigger. GRB 090709A 
turns out to have a rich phenomenology, featuring - in addition to 
the intrinsic photoeletric absorption - a clear spectral evolution (the 
photon index F steepens from 1.7 to 2.0 in ~2 days), as well as a 
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break in the flux decay, with a ~ 0.3 variation in the power law 
decay index occurring ~0.26 days after the burst. 

We will focus here on the extra absorption only, the evidence 
for which is fully confirmed - and strenghtened - by Swift data. In- 
deed, adopting a simple, redshifted neutral absorption model, all 
datasets yield consistent values for both the redshift and the red- 
shifted absorbing column, as well as of the spectral steepening as 
a function of time. XMM-Newton and Swift/XKI data point to a 
rather high redshift (z ~ 4.2) for GRB 090709A as well as to a 
huge column density in the host galaxy (JVh ~ 6 x 10 22 cm -2 ). 
Although a word of caution is required in considering such results 
(see Sect. 12.2.21 ), to our knowledge, this is the largest intrinsic ab- 
sorption ever observed in an X-ray afterglow spectrum (s ee e.g. 
Ide Luca et al]|2005l ; ICampana et alfcoOrj i lGrupe et~ai1l2007l) . 

This high absorption is also on the high side among expected 
column densities for sources located inside molecular clouds 
(Reichart & Price 2002). We note that part of the absorption could 
be due to intervening systems along the line of sight rather than 
being local to the host galaxy, as seen in the optical band in 
QSO and GRB studies ( the so-called Damp ed Lyman-a Absorbers; 
iReichart & Price] 1 20021 ; [Wolfe et all 120051) . Lack of optical spec- 
troscopy in the case of GRB 090709A prevents from drawing firm 
conclusions. 

Such a picture suggests an explanation of the peculiar, very red 
colors of the infrare d transient likely associated to GRB 090709A 
jMorgan et alj2009l) as due to reddening in the host galaxy. Indeed, 
if the GRB spectrum has no breaks in the X-ray to optical range and 
the redshift and intrinsic column derived from X-ray spectroscopy 
are correct, an extinction Ay ~ 3 in the GRB host galaxy (at 
z = 4.2) would fit the infrared data. This would point to a dust- 
to-gas ratio ~10% of that found in the Milky Way (assuming So- 
lar abundances), similar to the findings of o ther investigations (e.g. 
Galama &Wiiersll200ll ; lHiorth et al. 2003). Such an interpretation 
of the infrared data is possibly supported by the time decay of the 
infrared transient. The fading between the PAIRITEL and the Sub- 
aru observations is consistent with the t law describing the 
X-ray afterglow decay before the break at To + 0.26 days, which 
suggests a common origin for the X-ray and infrared emission. 

In any case, X-ray spectroscopy yields a robust indication that 
GRB 090709A was a very distant event. 

4.2 On the prompt temporal variability 

Coming back to the temporal properties of the prompt emission, our 
reanalysis, based on both Swift/BAT and INTEGRAL SPI/ACS data, 
could not confirm the presence of any significant signal at ~8.1 s. 
In any case, a complex multi-peaked light curve is clearly appar- 
ent, with a peak-to-peak separation of order ~8-10 s, and a peak at 
0.125 Hz is present in the power spectrum derived from BAT data, 
just below the 3a detection threshold. It is premature to comment 
on the apparent d iscrepancy between ou r results and those reported 
in a circular by Markwardt et al. ( 2009), in view of the lack of de- 
tailed information concerning the analysis performed by these au- 
thors. Indeed, several issues affect estimation of the significance 
of a signal when dealing with de-trending algorithms and Fourier 
transform techniques. Among the most important, we remember 
the presence of non-Poissonian noise in the power spectrum and 
the change in the statistical properties of a time series induced by 
operations such as subtraction and division. The former, if not taken 
into account properly, may result in an overestimation of the statis- 
tical significance of any peak sitting o n a non-Poissonian u nder- 
lying power spectrum continuum (see llsrael & S tella 1996). The 



latter affects the statistical properties of the original time-series. In 
particular, the division of a Possonian variable (the time series) with 
a (model-dependent) de-trending algorithm might result in a non- 
Poissonian distribuited time series with effects on the (unknown) 
statistical properties of the power spectrum noise. 

The claimed period of about 8 s is precisely in the range of 
periods of soft gamma-ray repeaters (SGRs) and anomalous X- 
ray pulsars (AXPs). These sources are thought to be magnetars: 
isolated ne utron stars powere d by strong (10 14 -10 15 G) magnetic 
fields (see [Mereghetti 2008 for a review). Thus, association of 
GRB 090709A with the activity of a magnetar was suggested, also 
possibly based on some similarity of the light curve to the one of a 
magnetar giant flare, featuring a sort of "pulsating tail". 

Based on high energy properties, we can rule out such an in- 
terpretation. First, we could find no unambiguous, coherent pulsa- 
tions in the 7-ray emission. Then, the spectrum of GRB 090709A is 
much harder than the typical spectra observed in the pulsating tails 
of the giant flares from SGRs. These are well described by thermal 
bremsstrahlung emission with kT ~ 15-30 keV, while the spec- 
trum of GRB090709A is typical for a GRB. Giant flares are also 
characterized by a short (<0.5 s), very bright and spectrally hard 
initial spike that was clearly absent in GRB 090709A (although this 
feature was lacking in some "intermediate flares" from SGRs, all 
of them had a very sharp initial rise, contrary to the light curve of 
GRB 090709 A). The location at high Galactic latitude (b ~ 20°) 
would also be very unusual for a Galactic SGR. Moreover, a Galac- 
tic origin may be safely ruled out by the observation of an absorbing 
column exceeding by a factor ~4 the Galactic value in the direc- 
tion of GRB 090709A. A giant flare from an extragalactic SGR can 
also be excluded, based on the absence of a visible nearby galaxy at 
the location of GRB 090709A. The observed fluence would imply 
a distance not larger than 600 kpc (assuming an energetic similar 
to the one of other observed magnetar giant flares). The enormous 
energy requirement implied by a cosmological distance rules out 
the SGR giant flare hypothesis. 

As already discussed in the previous section, we conclude that 
GRB 090709 A is a standard, long GRB, with a multi-peak struc- 
ture in the prompt emission. We will not go into pure speculations 
by considering physical processes that could produce a true quasi- 
periodic signal. 

The variable prompt emission could be the signature of non- 
stationary processes in the GRB in ner en gine. As a likely possi- 
bility - discussed by Bes kin et al.l d2009l) who observed a simi- 
lar phenomenology in the optical emission from the "naked-eye" 
GRB 0803 19B (but with a possibly stronger evidence for quasi- 
periodicity) - the peculiar variability could be related to cyclic 
accretion by the central newborn compact object. Such a phe- 
nomenon could be due to the fragmentation of an accretion disc 
due to some kind of instability. For instance, iMasada et al.l d2007h 
explain short-time variability in the prompt emission by GRBs as 
due to magneto-rotational instabilities developing in a massive, hot 
hyperaccreting disc surrounding a central black hole of a few stel- 
lar masses. At a redshift z ~ 4.2, the variability time scale in 
th e source frame woul d be of ~ 1.5 s. Thus, adopting the model 
of IMasada etafl j2007h , the observed properties of GRB 090709A 
could fit into such a scenario by assuming a beaming factor ~100, 
for a ~1 Mq accretion disc with inner radius of ~30 gravitational 
radii, surrounding a ~4M© central black hole. 
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